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We study 1-D diffusion of TV hard-core interacting Brownian particles driven by the space- and 
time-dependent external force. We give the exact solution of the A'^-particle Smoluchowski diffusion 
. . , equation. In particular, we investigate the nonequilibrium energetics of two interacting particles 

Cn ' under the time-periodic driving. The hard-core interaction induces entropic repulsion which differ- 

entiates the energetics of the two particles. We present exact time-asymptotic results which describe 
^^ , the mean energy, the accepted work and heat, and the entropy production of interacting particles 

^^ ' and we contrast these quantities against the corresponding ones for the non-interacting particles. 

^- 

^ . PACS numbers: 05.40.-a, 87.10.Mn, 87.15.hg, 02.50.Ga, 05.40.Jc, 05.10.Gg, 82.75.Jn 
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Introduction: Stochastic dynamics of interacting particles in a one-dimensional environment is both of great practical 
and theoretical interest. Due to the one-dimensionality of the problem, the inter-particle interactions play a crucial 
role as they alter qualitative features of the particle dynamics. The type of interaction we deal with in this Letter is 
a so called hard-core interaction. 

The diffusion of particles in narrow channels, where the particles cannot pass each other and their relative ordering 
is conserved, is known as the single-file diffusion (SFD). The concept of SFD was first introduced by Hodgkin and 
Keynes in relation to the transport of water and ions through the molecular-sized channels in membranes 1]. Since 
then, numerous examples of SFD in biological, chemical, and physical processes were studied (e.g. transport of 
adsorbate molecules through zeolites with a one-dimensional channel system 0, Q, geometrically constrained nano- 
sized particles in nano-sized pores [Jl , migration of adsorbed molecules on surfaces 0j , diffusion in nanotubes [yJH, a 
carrier migration in polymers and superionic conductors [8| , diffusion of colloids in one-dimensional channels |9l4lll| , 
confined dynamics of millimetric steel balls [l_2|). 

While the global properties of a SFD-system are identical to those for the system of independent particles, the 
dynamics of an individual particle (also called tagged particle or tracer) is considerably different [l3l - [l5l | . Theoretical 
description of SFD was introduced by Harris in 1965. In his pioneering study [16| he showed that the mean square 
displacement of the tagged particle increases with time as i^ " (in contrast to its linear increase for the single free 
particle) . This result was subsequently reestablished by many other authors using different mathematical tools (for a 
comprehensive review cf. the Introduction in |17|). The first exact solutions of the diffusion equation for SFD systems 
appeared only recently. The solution for an arbitrary number N of identical particles diffusing along the infinite line 
has been obtained by Rodenbeck et al. in [Ij] via the reflection principle. Using a different theoretical procedure, the 
result has been re-derived in [18|. The exact solution for the diffusion within the finite 1-D interval has been found 
in [13] through the Bethe ansatz. The exact solution for N = 2 particles having different diffusion constants has 
been obtained in [13|. Another remarkable exact result for A^ = 2 has been derived in [2^ where the authors assume 
the finite-range interaction between the particles. The first calculation including the time-independent external force 
has been published in [2l[ . Using a so called " Jepsen line" , the probability density for the single tagged particle was 
obtained in the limit of an infinite system with the constant particle density. As far as we know, no exact solution of 
the SFD of A^ particles driven by the external time-dependent force has been published, yet. 

The main objective of this Letter is to analyse the thermodynamic properties of an exactly solvable, experimentally 
verifiable model of an externally driven SFD. Hence the setting must incorporate the following attributes: a) the 
hard-core interaction among the particles, b) the non-trivial equilibrium state if the external driving is switched off, 
c) the time-periodic external driving inducing a periodic response. We first give the exact solution ^ of the diffusion 
equation with an arbitrary time- and space-dependent external driving of the A^ hard-core interacting particles. 
The solution emerges in the form of a sum of products of the single-particle probability densities which solve the 
corresponding one-particle diffusion problem. Differently speaking, in the second step, we need an exactly solvable 
one-particle model which would include the above ingredients b) and c). We focus on the problem which is related 
with the famous barometric formula. We consider a particle diffusing in a half-space x E (0, +oo) and acted upon by 
the spatially homogeneous and time-dependent force with a refiective boundary placed at x = 0. Using the general 
A^-particle solution and the single-particle probability density for the barometric problem we investigate in detail the 
diffusion of two interacting particles. 

General solution: Consider A^ identical hard-core interacting Brownian particles [22l | , each with the same diffusion 
constant D, diffusing under the influence of the external force F(x, t). Let the vector x = xi, ... , xn denotes positions 
of the particles at the time t and let y — yi,... ,yN stands for their initial positions at to < t. The time evolution 
of the joint conditional probability density for the positions of the particles, say p^^\x] t \ y; to), is controlled by the 
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with the initial condition 
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In Eq. ([T]) m and 7 are the mass of and the friction coefficient for the single particle, respectively. Suppose that the 
particles are initially ordered as j/i < 2/2 < ••• < yN- Owing to their hard-core interaction, the ordering is conserved 
at all t > to and the function p^^'{x;t\y;to) vanishes outside the domain a;i < 2:2 < ... < xn- This restraint is 
guaranteed by the non- crossing conditions 
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j — 1, ... ,N — 1. In order to derive the conditions, one first introduces the A^-particle probabihty current J — 
(Jf , ..., j]^ ) with the j-th component 

jf)(f;i|y;M- I^^F{x„t) ~ D^^ p^''\x;t\yM . (4) 

Thereupon the hard-core interaction can be implemented through the requirements 

jj^' (f ; 1 1 y; to) ~ jj^i' {x; t \ y; to) = . (5) 

Xj-^l — Xj 

Inserting the components (|4]), we arrive at the conditions ([3]). Notice that the argumentation holds true only for the 
identical particles. In other words, the parameters D, m, 7, and the external force F(x,t) must be the same for any 
individual particle from the SFD system. 

Let us now consider the Smoluchowski diffusion equation ([Ij with A^ = 1. Assume we know the solution 
p^^'{x;t\y;to) of this single-particle diffusion problem. Then we claim the following. The exact TV-particle joint 
probability density p'^^\x]t \ y;to) which fulfils (HI), ^ and ^ reads 

N\ 

p(^'(f; 1 1 y; to) = ^p'^'^\xi;t \ y-„^(i);to)p^'^\x2;t\ ynki2);to)--p''^\xN;t \ y^,,(Ar); to) , (6) 

fc=i 

if xi < a;2 < ... < xat, and it vanishes otherwise. The summation is taken over all permutations TTfc. A given 
permutation specifies the initial conditions of the single-particle densities on the right hand side. The claim can be 
proven by a direct check: we simply insert the right hand side of Eq. © into ^ and ([3]). 

The formula ([6]) expresses the exact solution of the many-particle problem with the hard-core interaction through a 
simpler object which is the single-particle probability density. Notice that in the special case of the SFD without the 
external driving {F{x, t) = 0) the result ^ correctly reproduces the solutions which were previously derived using 
the refiection principle [1J|, or the Bethe ansatz 17]. However, these techniques fail in the presence of the external 
time-dependent force F{x, t). This is the reason why the SFD under the external driving has remained unsolved, yet. 
Our derivation demonstrates that the structure of the solution as given by Eq. ^ remains valid even in the SFD 
under the external driving. 

Barometric SFD of two particles: Place two Brownian particles into the long 1-D channel, plug one end of the 
channel and let each particle be acted upon by the space- homogeneous and time-dependent force F(t). In other 
words, let the two particles diffuse in the external potential 

{- Xi F{t) for Xi> , 
(7) 
+00 for Xi < , i — 1,2 . 

The driving force incorporates two components, 

F{t) ^ Fo + Fi sin(wi) . (8) 

The time-independent component Fo alone would push the particles to the left towards the reflective boundary (if 
Fq < 0), or to the right (if Fq > 0). The time-dependent component Fism{ujt) harmonically oscillates with the 
angular frequency lo. Having in mind this scenario, we wish to contrast the dynamics and the energetics of the system 
of two interacting particles against the model without the interaction. In the latter case the analysis trivially follows 
from the solution concerning the single- diffusing particle. 

Assuming the driving force (|S]), the most interesting physics emerges if the oscillating component Fi sin(a;i) super- 
poses with a negative static force Fq < 0. In this case, Fq acts against the general spreading tendency stemming from 
the thermal fluctuations of the surroundings. Due to the periodic driving the system of particles approaches a definite 
steady state exhibiting cyclic changes of the probability density regardless the initial positions yi , j/2 ■ A necessary 
and sufficient condition for the existence of such steady state is Fq < [24]. Let us denote probability density of 
the single-diffusing particle in the steady state as ps{x;t) (i.e., in the steady state, p^^'{x;t\y,to) ~ ps{x;t), and 
ps{x; t -\- 27r/w) = ps{x; t), for any x and t). In the steady state, the right hand side of the formula (O collapses into 
the product of just N terms. Explicitly, for N — 2, the two-particle joint probability density in the steady state reads 

p^'^\xi,X2;t) = 2ps{xi ; t)ps{x2 ; t) , (9) 



if Xi < X2, and it vanishes otherwise. Hence the marginal probability densities for the position of the left (L) and the 
right (R) interacting particle are 
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respectively. Carrying out the integrations and using the normalization condition for the function ps(x;t), we have 
proven an important identity 

PR{x;t)+pL{x;t) = 2ps{x;t) . (12) 

For any t, there exist a unique coordinate £,{t) such that pl(^(^); t) = PR{£,{t)', t)- For x ^ £,{t), Pr{x; t) ^ Pl{x; t). 
Being a periodic function of time, the probability density psix] i) can be represented as the Fourier series 
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The coefficients Uk{x) were found in '24'|. Using the standard matrix notation, they assume the form 

Uk{x)=mj^-^{k\LE{x)R+\f) , 



(13) 
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where |/) is the column vector of the complex amplitudes, fk ~ (k |M_^ | O), and the matrices E(a::), L, M± possess 
the matrix elements |23| 
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where to and n are integers, k, ~ |Fo|Fi/[2aj_D(TO7)^], C — 4:UjD{m'^)'^/FQ, and Ife(-) stands for the modified Bessel 
function of the first kind, of the order k. Finally, introducing the coefficients 
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the marginal densities ([TU]), pT|) can be written in the forms 
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These formulae, together with Eq. ([131) . represent an exact asymptotic result and they form the basis of the further 
discussion. In the numerical illustrations (Fig. 1) we had to reduce the infinite matrices (|15l) . (IT51) . p7)) into their 
(finite) central blocks. 

Mean position represents the centre of mass of the probability density. In the steady state, we define mean positions 
of the individual particles as 
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Due to the oscillatory driving ([5]), the mean positions will oscillate with the fundamental frequency w. In an average 
sense, the hard-core interaction produces a repulsive force among the particles. It (in average) shifts the right (left) 
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FIG. 1. The driving force ((8| (panel (a)); mean works 1)25^ (panel (b)); mean energies (|23|l (panel (c)); and mean heats (|27p 
(panel (d)) as the functions of time in the steady state. In (b), (c) and (d), the solid black line, the blue dashed line, and the 
red dot-dashed line correspond to the single-diffusing particle, the right particle, and the left particle, respectively. We have 
usedFo = -l.ON, Fi = 3.0N, a; = 0.4 7rs~\ D = 1.0m^s~\ m,-) = 1.0kgs~\ 



particle to the right (left) as compared to the case without the interaction, i.e., /XL(i) < Ms(^) < Mr(0 holds at any 
instant t. 

If we average the external potential (/)(x, t) over all possible positions of the particular particle at a given instant 
we obtain the mean energy of that particle: 



Ea (t) = -[Fo + Fi iim{ojt)]fia {t) , a == S, L, R , 



(23) 



Generally speaking, the internal energies Es{t), Ei^{t), E-[i{t) are periodic functions of time with the fundamental 
period 27r/w. Their oscillations express the combine effect of both the periodically modulated heat flow to the 
surroundings and the periodic exchange of work done on the particle by an external agent. From Eq. (IT^ . it follows 
that the total mean energy of two interacting particles is equal to the total mean energy of two non-interacting 
particles, i.e.. 



Ei^it) + En{t) ^ 2Esit) 



(24) 



In other words, the hard-core interaction does not contribute to the total energy. Therefore, the repulsive force among 
the particles arises from the purely entropic effect. On the other hand, the mean energies of individual interacting 
particles are changed significantly as compared to the case without interaction (see Fig. [T]). 

The mean work done on the particle by an external agent during the time interval [0, t] reads [25| 



Wait) = -Ficj dt' cos{ujt')fia{t') , a = S,L,R. 
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The total work done on the system of two interacting particles again equals to the total work done on two non- 
interacting particles. 



WLit) + WRit)^2Ws{t) . 



(26) 



Nevertheless, as was stressed within the discussion of mean positions, the hard-core interaction in average shifts the 
right (left) tagged particle to the right (left) as compared to the case without interaction. Hence the absolute value 
of work done on the right (left) particle is always bigger (smaller) than in the case without interaction. Since in our 



setting the diffusion is just a non-equilibrium isothermal process, the total work done on the particles per one period 
27r/w is always positive. This work is entirely dissipated into heat. However, if we allow for a temperature modulation 
and choose an appropriate temperature schedule, the system can act as a heat engine which converts a part of the 
heat accepted form the environment into the useful work. If this is the case, the right particle will perform a bigger 
work as compared both to the left particle and to the single-diffusing particle. 

The mean heat released to the environment and the total entropy increase per period are intimately related. From 
the first law of thermodynamics it follows that 

Qcit) = -[E^it) -Wc.it)] , a = S,L,R. (27) 

is the heat dissipated during the time interval [0, t] by the individual particles. The entropy generated in the environ- 
ment due to the dissipative motion of the individual particle during the time interval [0,i] is then Sa{t) = Qa{t)/T. 
Having periodic changes of the internal energy, the heat dissipated during one period is equal to the work done on the 
system during one period Qa{2Tr/uj) = Wa{2Tr/uj). Divided by the temperature of the surroundings we get the total 
entropy increase per one period Sai^n/uj) = Qai2iT/uj)/T, and again S'L(27r/a;) -|- 5'R(27r/a;) = 25's(27r/w). That is, 
the entropic repulsion stemming from the hard-core interaction does not influence the total entropy production. 

In summary, the global (aggregative) quantities (e.g. the total mean energy (24), the total mean work (26), the 
total mean entropy production per period) evaluated for the system of interacting particles equal to those for the 
system of independent particles. Physically, the conclusion derives from the zero range of the interaction among the 
particles. The observation holds for a general external driving and for an arbitrary number of particles, N. Contrary 
to this, the simple contact interaction strongly influences the one particle dynamical characteristics. 
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